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ABSTRACT 

Structure formation provides a strong test of any cosmic acceleration model because a suc- 
cessful dark energy model must not inhibit or overpredict the development of observed large-scale 
structures. Traditional approaches to studies of structure formation in the presence of dark energy 
or a modified gravity implement a modified Press-Schechter formalism, which relates the linear 
overdensities to the abundance of dark matter haloes at the same time. We critically examine the 
universality of the Press-Schechter formalism for different cosmologies, and show that the halo 
abundance is best correlated with spherical linear overdensity at 94% of collapse (or observation) 
time. We then extend this argument to ellipsoidal collapse (which decreases the fractional time of 
best correlation for small haloes), and show that our results agree with deviations from modified 
Press-Schechter formalism seen in simulated mass functions. This provides a novel universal pre- 
scription to measure linear density evolution, based on current and future observations of cluster 
(or dark matter) halo mass function. In particular, even observations of cluster abundance in a 
single epoch will constrain the entire history of linear growth of cosmological of perturbations. 

Subject headings: cosmology: large-scale structure of universe 



1. Introduction 

An important part of the effort to explain cos- 
mic acceleration and the cosmological constant 
problem is testing proposed models in the context 
of what are, at this stage, better-established phys- 
ical pictures. Structure formation could prove to 
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be an incredibly useful phenomenological method 
for distinguishing models of cosmic acceleration. 

It is currently believed that large-scale struc- 
ture formation has its seeds in sma ll quantum fluc- 
tuatio ns in the early universe fe.g-. lMukhanov et al 
(|1992l )). The current model for structure forma- 
tion is elegant in its fundamental simplicity. Ran- 
dom inhomogencitics, artifacts of cosmic inflation, 
create a runaway effect in which overdense regions 
attract more matter, thus becoming more dense 
and leading to galaxies, stars, and planets. 



1 



Better understanding of this process is inde- 
pendently an intriguing enterprize in the field of 
cosmology. However, in this work we focus on 
the relationship between the cosmic acceleration 
and structure formation. More specifically, differ- 
ent cosmological pictures (cosmologies with differ- 
ing causes of acceleration, such as a cosmologi- 
cal constant, dark energy, or modifications of Ein- 
stein gravity) might have expansion histories that 
are similar to one another but leave different im- 
prints on large-scale structures, and in particular 
on galaxy clusters. Therefore, structure forma- 
tion provides a unique testin g ground for mod - 
els of cosmic acce l eratio n fe.g.Jlshak et al ] (l2006l) : 
Acauaviva et all (|2008l ): IVikhlinin et al.l (|2009t )i. 
Here we critically examine some of the assump- 
tions in this program, and develop a framework to 
enhance the accuracy of this kind of work. 

The first step in this directi on is to revisit how 
the Prcss-Schcchter formalism (Press fe Schechter 



( 20081 ) propose a fitting formula of the form: 



1974 ) (PS F) is used to pred i ct the cluster mass 
function. IPress fc Schechter ( 1974 ) have argued 
that the number density of dark matter haloes (or 
galaxy clusters) of mass M is given by: 



dn(M, z) 
dM 
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where a 2 (M,z) is the variance of linear overden- 
sity in spherical regions of mass M at redshift z, 
while p m (z) is the mean matter density of the uni- 
verse. For random gaussian initial conditions, f[a] 
is given by: 
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where 6 SC (~ 1.68 for most ACDM cosmologies) is 
the sp herical collapse thre shold for linear overden- 
sities (jQunn fc QottJI 19721) . 



While the PSF successfully predicts the broad 
features of the simulated cluster mass functions, 
it proves too simplistic for detailed model com- 
parisons required for precision c osmology. Conse- 
quent l y, several authors incl uding Sheth fe Tormen 



Evrard et al l (120021) . 
Tinker et all (|2008l ) 



( 19991). I Jenkins et all (120011). 
Warren et al.l (|2006f ). and 
have refined the function / (er) to better match the 
mass functions in N-body simulations. We shall 
refer to t hese approach e s as modifi ed PSF. For 
example, IWarren et al.l ([20061 ) and iTinker et al 
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where (A,a,b,c) = (0.186,1.47,2.57,1.19) give a 
good fit to simulated haloes of overdensity A = 
200, in a concordan ce ACDM cosmology at z = 
([Tinker et alJl2008t ) (see Fig. [T}. While most of 
this work is bas ed on fitting f ormul ae to simulated 
mass functions, ISheth et al ] ([20011 ) argue that an 
approximate implementation of ellipsoidal collapse 
can account for most of the deviations from the 
PSF. 
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Fig. 1. — A comparison of Prcss-Schechter prediction for 
the function /(c) (dotted line; Eq. with a parameterized 
fit to the numerical simulations (solid line; Eq. [3J. 

However, a more pressing question for cosmo- 
logical applications is whether the function /(a) 
is universal, or rather can vary for different cos- 
mologies or cosmic acceleration models. In other 
words, could the same modified PSF be used to 
accurately predict halo abundance in cosmologies 
with different cosmological parameters? While 
earlier studies failed to find such dependence, 



Tinker et al.l (J2008) first noticed a systematic evo- 
lution of /(cr) with redshift, implying a brcak- 
down of universality at the 1 — 30% level (also see 
Bhattacharva et al l (|2010h ). ICourtin et all (|2010l ) 
note that universality is limited by the nonlinear- 
ity of structure formation, and the cluster mass 
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function shows some redshift dependence at higher 
redshifts that can be (partially) understood in the 
context of spherical collapse. However, the spher- 
ical collapse model falls short of explaining the 
10 — 30% deviations from universality in all but 
the most massive clusters (see Fig. |6j). 

In this work, we contribute to the effort to bet- 
ter understand the role and limits of universal- 
ity in the cluster mass function by introducing a 
new parameter that appears to be universal across 
cosmological modelsLj In particular, the modi- 
fied PSF relies on c(M), the root- mean-square 
of linear density fluctuations at the time of ob- 
servation, when in reality, observed clusters are 
very non-linear objects with overdensities exceed- 
ing 200. We thus seek to find a universal time 
in the past when we could make a connection be- 
tween the nonlinear structures that we observe in 
the present and the linear structures that existed 
in the past, since all structures go through a linear 
phase. Our basic strategy is to find the time in the 
past at which the linear density of the structures 
that collapse today show minimum dispersion, as 
we vary cosmologies. 

In § [21 we introduce linear perturbation theory 
for modeling structure formation which leads to a 
linear differential equation. 

In § El we present the complete nonlinear differ- 
ential equation that governs the growth of matter 
perturbations in spherical overdense regions in the 
presence of a cosmological constant. 

In § 21 we describe a numerical method which 
we developed in order to solve both the nonlin- 
ear and linear structure formation equations in 
the presence of a cosmological constant, or dark 
energy with a constant equation of state. 

In § we discuss the implications of our nu- 
merical study, i.e. we find a universal fraction of 
collapse time, where linear density is the same in- 
dependent of cosmology. 

In § [6] we reexamine the idea of universality of 
mass functions in light of the results of the previ- 
ous section, including the effect of ellipsoidal col- 



1 Here we define cosmological models to mean different val- 
ues of the cosmological constant density Q\ or different 
values of the dark energy equation of state parameter, w. 
For now, we do not consider quintessence models where w 
is dynamical in nature, or possible modifications of Ein- 
stein gravity. Future work will extend our study to such 
possibilities. 



lapse on the formalism, and propose a new mass 
function for general dark energy models. 

Finally, in § [71 we conclude with an overview of 
our results and a discussion of future prospects. 

2. Background: Linear Perturbations 

The linear perturbation theory that is used to 
describe structure formation can be derived via 
a fluid picture. We use a Newtonian treatment 
because when density perturbations are small, 
the gravitation al potential will be nonrelativistic 
([Peebles! 119931 ). The standard equations of fluid 
dynamics are reviewed here. First, the continuity 
equation: 



| + v.(^) = o ) 



(4) 



where p is the matter density and v is the fluid 
velocity. The Euler equation is: 

dv ^ „ w VP „ J 

¥ +«'V« = V0, 5 

at p 

where 4> is the gravitational potential. Finally, the 
Poisson equation is: 



V 2 </> = AnG N p. 



(0) 



Introducing the comoving coordinates x = aCR' 
we can write the velocity in terms of the comoving 
coordinates and the scale factor a(t): 



dr 
dt 



-r + u(-,t) (7) 
a a 



where u = xa(t) is the peculiar velocity. Note that 
at constant r, 4| = so ax + ax = 0, giving us 
^x = — ^r. This leads us to the following relation: 
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dt 
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We now wish to rewrite Equation 2] in terms of 
the comoving coordinates, which essentially means 
replacing the partial differential with the modified 
one from Equation [51 Moreover, what we are re- 
ally interested in is the development of relative 
deviations from the mean density, or the density 
contrast: 5 = = — 1. Thus, we write p = p(l + 5) 
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and assuming that p is the mean density of regular 
matter, we expect p oc a~ 3 . This gives: 







+ iv.(l + *)* (9) 



We make similar transformations for the Poisson 
and Euler equations. Next, we drop higher order 
terms (e.g. 0(u 2 ) or uS). We differentiate the 
linearized continuity equation and take the diver- 
gence of the linearized Euler equation. This gives 
us a differential equation that depends entirely on 
5 and not on u: 



2 a -5 
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(10) 



The linear growth factor D(t) is defined as the 
growing solution for 5 in this equation. 

3. A &: Non-linear Structure Formation 

In § [21 wc derived the differential equation that 
governs the growth of linear matter perturbations. 
We used a familiar fluid dynamics picture to do 
so. Here we derive the full non-linear equation 
for spherical overdensities using only cosmological 
considerations. It should be noted that this partic- 
ular form of the non-linear equation is only strictly 
valid for ACDM cosmologies, where the inside of 
a spherical top-hat overdensity can be considered 
as a separate universe. More complex models such 
as dark energy models with different values of w 
require additional considerations which will be dis- 
cussed below. 

We consider a physical picture in which struc- 
ture formation arises due to a uniformly positive 
spherical perturbation away from an average mat- 
ter density, i.e. a top-hat matter overdensity. This 
scenario is similar to considering two cosmologies 
with two distinct scale factors: one for the outer 
universe and another for the overdense region. Of 
course, we are interested in a scenario where a 
dark energy component similar to a cosmological 
constant is at play, so we will assume the presence 
of one as part of our base model. 

For the external universe, we write the Fried- 
mann equation with zero curvature: 
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As we did in § [2l we can assume p oc a~ 3 for or- 
dinary matter, while p\ = const, denotes the cos- 
mological constant density. We also note that the 
value of the cosmological constant will be the same 
inside the overdense region and the background. 

In a general scenario with dynamical dark en- 
ergy models we cannot assume that curvature, of- 
ten denoted by k, will be a constant inside the 
overdense region due to the presence of p r essure 
gradients. Therefore, IWang fc Steinhardt (Il998h 
point out that we are compelled to use the time- 
time component of Einstein's equations, as these 
do not explicitly involve the curvature term. How- 
ever, in the presence of a cosmological constant, or 
w = — 1, we can ignore these considerations and 
begin with the first order Fricdmann equation. We 
then calculate k, which can be seen as an integra- 
tion constant that arises in going from second to 
first order Friedmann formulations, using initial 
conditions. 

The scale factor in the overdense region is gov- 
erned by: 
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Again, we define pi = p(l + 5). A little bit of al- 
gebra yields the following full differential equation 
for 8: 
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(13) 

It is important to reiterate the importance of 
having access to both linear and n onlinear solu- 
tions. As noted bv lPace et all (|2010h amongst oth- 
ers, although initially both the complete solution 
and its linear approximation will track, eventually 
the nonlinear solution will grow much faster rela- 
tive to the scale factor. 

Following iLiddle fc Lvthl (|2000h . we can show 



-k = 



that knowing nonlinear theory is necessary. A crit- 
ical point in the evolution of a structure's collapse 
is the turnaround event in which universal expan- 
sion's dominance over the perturbation is eclipsed 
by gravitational collapse. In other words, at the 
turnaround point, a potential structure has de- 
tached from background expansion, but complete 
gravitationally-bound structure formation has not 
yet begun. This might be thought as the true birth 
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of a structure within the void. Knowing the non- 
linear solution allows us to find out the value of the 
scale factor and the overdensity at this so-called 
"turn-around point." 

Numerically, at the point of complete collapse 
the nonlinear solution will "blow up" and ap- 
proach infinity (see Fig. [5]). The linear density at 
this point, S sc is the quantity that enters the origi- 
nal Press-Schechter formalism (see §[T]) and is used 
as proxy between linear and non-linear structures. 
Physically, we do not expect an actual singularity. 
This "blow up" point is considered to be the be- 
ginning of virialization, a process whereby energy 
in the bulk infall of matter is redistributed into 
random motion of dark matter particles, leading 
to a system in virial equilibrium, where kinetic en- 
ergy T and the potential energy U are related by 
the virial theorem: 



T 



1 



(R 8U/dR) x 



(14) 



see aor fc Lahavi (|2005D ). For the purposes of 
this study, more details on this process are not 
necessary. 

We can generalize to cases where w ^ — 1, i.e. 



non-co smological constant scenarios. lAbramo et al 
(J2007|) provides a comprehensive derivation of the 
full non-linear equation, which we refer the in- 
terested reader to for complete details. For our 
purposes, it will suffic e to show the final r esult, 
which is Equation 7 in lAbramo et al] (|2007j) : 
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l + WjJ' J 3(1 + ^)1 + ^ 



-AttGY, PkO- + «>k)(l + 3w k )5 k (l + 4) = 0. (15) 

k 

The subscripts j and k refer to different fluids in 
the system, e.g. matter and dark energy. In a 
scenario where the dark energy can clump, this 
becomes a system of equations. As we discuss be- 
low, we do not allow this possibility. Therefore, 
noting that w = for matter, we get the follow- 
ing non-linear equation governing the behavior of 
matter perturbations: 

4 Si 



3 1 + <L 



- 47rGp m S m (l + S m ) = 0. 

(16) 

But, why exactly are we allowed to ignore the 
clumping in dark energy? 



The discussion about dark energy perturbations 
is often cast in terms o f the effective speed o f sound 
for dark energy, e.g., Bean fc Pore ( 20041 ). Typ- 
ically one may expect that for an adiabatic fluid 
with a constant equation state parameter w, the 
speed of sound is given by c 2 = SP/5p = w. How- 
ever, when w c ie < 0, clearly this becomes imagi- 
nary, suggestin g a cat astrophic instability, e.g., see 



Afshordi et all (|2005h . and thus we must use a 
more general definition of the c. In other words, 
dark energy with constant w cannot be an adi- 
abatic fluid, implying w = — ^ j 2 . While the 
equation of state parameter can still give us infor- 
mation about the background evolution, the full 
action of the fluid is necessary to compute its effec- 



eff 



|2. It turns out that 

Sp 



tive speed of sound: c; 

for the simplest quintessence models c e // = 1, al- 
though for more general actions c e ff could essen- 



Armendariz-Picon et al 



tially take any value, e.g 
(|200lh . 

As to the question of dark energy clumping, 
we know that pressure fluctuations propagate with 
the speed of sound c e //. Therefore, dark energy 
should be smooth on scales smaller than its sound 
horizon ~ c e ff/H. As long as c e // ~ 1, all the 
collapsed structures at late times are much smaller 
than the sound horizon, implying that dark energy 
perturbations 5de should be negligible for their 
formation. 

4. Numerical Techniques 



Although iGunn fc Gottl (|1972l ) showed that 
the perturbation equation with the assumption 
of spherical symmetry (also known as the spheri- 
cal top-hat problem) can be solved analytically for 
the case of an Einstein-de Sitter (EdS) universe 
(A = 0), we are interested in cases where the 
cosmological constant/dark energy are non-zero. 
Therefore, a numerical solution is necessary. 

We built a code in C++ that utilizes the Rungc- 
Kutta method for numerical solutions of ordinary 
dif ferential equations. The f ull code can be found 



Prescod-Weinstein <l201dh . The code runs an 



instance of a loop for value of 51a 



present 



in the 



range of to 0.7 (the currently measured value 
of dark energy density), which solves the second- 
order linear differential equation, Equation I1Q[ as 
well as the full non-linear equation, Equation 1131 
In order to solve both equations, the solution to 
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the differential equation for the background scale 
factor, Equation II 1| is found. The curvature con- 
stant is calculated for the initial overdensity be- 
fore proceeding to a solution for both the linear 
and nonlinear differential equations. 

Collapse time is formally defined to be the time 
at which nonlinear overdensity, > goes to infin- 
ity. However, we define collapse numerically by 
requiring a large value of nonlinear overdensity, 
Snl — 200. The initial conditions are calibrated 
such that they provide the same results as the ana- 
lytical EdS model, namely that SLitcollapse) ~ 1-7. 
This is essentially done by assuming that at early 
times 8 scales linearly with the scale factor a, as 
expected for the growing solution to Equation ([TU)) 
in the matter era^. 
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Fig. 2. — This plot shows nonlinear vs. linear overdensi- 
ties in the Einstein-de Sitter Universe. The code produces, 
as expected, a &nl that diverges when 8l ~ 1.68. 



5. Results 

As stated in the Introduction, we wish to dis- 
cover the time (as a fraction of the collapse time) 
such that the variance of <5l(£) is at a minimum, 
as we vary cosmologies. In doing so, we find the 
time in cosmological history where linear theory 
is most likely to accurately predict gravitational 
collapse, independent of cosmology. 

For simplicity, we chose to plot the relative dif- 
ference of 5 Lit) for different cosmologies with re- 
spect to the fraction of time to the collapse time 
in the Einstein-de Sitter universe. However, it 
should be noted that the results arc independent 
of this choice, and one could easily calibrate with 



respect to a universe with a non-zero A instead. 
Having found a common ground, all the data was 
searched for a single point in time (in units of col- 
lapse time) where the variance of g^jgL — 1 was 
at a minimum, where Q stands for different dark 
energy models under consideration, whether a gen- 
eral fluid or a cosmological constant. 

The interpolations and variance computations 
were carried ou t using a script in MatLab , which 
can be found in Prescod-Weinstein ( 2010l) . In the 
case of the simple cosmological constant models, 

- 0.94, the 



taolla 



Fig. [3] shows the result that at 
variance in SL{t/t co iiapse) for different cosmologi- 
cal constant models hits a minimum of 1.80 x 10~ 9 . 
Values of 6l at range from 1.602 for Einstein-de 
Sitter to 1.599 in the A = 0.7 universe. 

Computing the solution to Equation (|16|) re- 
quires some modifications to the code. Instead of 
scanning over different cosmological constant val- 
ues, this version of the code varies between con- 
stant values of w. Moreover, as current observa- 



tions, e.g.. iKomatsu et ail ( 2010| ) set a (very con- 
servative) upper limit of —1/3 for the value of w , 
we studied cases where w was smaller than —1/3. 

Because the cosmological constant is a special 
case of this scenario with w = —1, we are able 
to check the self-consistency of our methods (Eq. 
IT51 vs. Eq. IT6"| finding that both versions of the 
code produce the same results for a universe with 
JIa = 0.7 (approximately the universe that we live 
in). 

Fig. (|4j) shows a similar comparison to that of 
Fig. ([3]), but with different values of equation of 
state, w, with VLde = 0.7. Interestingly, we can 
again clearly see a point of minimum variance at 
t/tcoiiapse — 0.94, suggesting that this result might 
be quite independent of the dark energy model (at 
least within the spherical collapse approximation) . 
The variance at this minimum is 7.5 x 10" 7 . In the 
next section, we discuss the physical significance 
of this result. 



3 This is below the predicted analytic value of 1.68 at collapse 
because while 8l = 1.68 is expected as &nl — > °o> we 
have set the collapse to occur at <5jvz, = 200, resulting in a 
lowered collapse 8^. 



2 It is important to remember that 8 is a ratio of two numbers 
whose units are that of density. Therefore, 8 is dimension- 
less. 
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Fig. 3. Relative change in Sj^ft / ^collapse 

) in ACDM 

cosmologies (f2y\ = 0.1, 0.2, 0.7), with respect to the 
Einstcin-de Sitter Universe. t co u apse is calculated for 
spherical overdensities. The curves seem to intersect at 
t/tcollapse — 0.94, and a calculation of the point of mini- 
mum variance between the lines confirms this. 




Fig. 4. — Same as Fig. (|3j) , but with varying values 
of dark energy equation of state (ui = —0.3, —0.4, —1), 
at Q^(today) = 0.7. Again, curves seem to intersect at 
collapse — 0.94, and, again, a calculation of the point 
of minimum variance between the lines confirms this. We 
begin at the value -0.3 because current data constrains the 
parameter to be smaller than this value. 



6. What does the Cluster Mass Function 
teach us about cosmology? 

At first, one might be puzzled by the fact that 
Sl happens to have almost the same value at 94% 
of the collapse time, independent of if or A, even 
though the linear approximation breaks down long 
before this point. In other words, why should non- 
linear collapse show such strong correlation with 
the linear evolution? This can be understood as 
the near cancelation of two different effects with 
opposite signs: 

1. With the exact same initial conditions, the 
presence of dark energy weakens the gravita- 
tional attraction near the turn around point, 
which in turn stretches the collapse time. 



2. The linear growth factor D(t), which is the 
growing solution to Equation (|10|) . slows 
down as dark energy starts to dominate, 
since the Hubble friction 2HS stops decaying 
(H — > const.), while matter density decays 
exponentially p cx a~ 3 oc exp(— Sift). 

It turns out that for near ACDM cosmologies, 
these two effects nearly cancel each other, i.e. the 
linear growth is slowed down, but t co u apse is also 
longer, resulting in nearly the same value of Sl 
close to (i.e. at 94% of) the collapse time. 
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Fig. 5. — The expected time of minimum variance in 
linear overdensity, in units of ellipsoidal collapse time. Ob- 
serving the cluster number counts at t co i should tell us the 
linear overdensity of the collapsing region at t», indepen- 
dent of the dark energy model. 

However, we should note that this result is 
specific to the spherical collapse scenario. On 
the other hand, real collapsing regions are far 
from spherical. Even though very rare peaks 
of a random gaussian field can be approximated 
as spheres, more abundant haloe s could signifi 



cantly deviate from sphericity, e.g.. lBardeen et al 



1986). For Einstein-de Sitter universe. ISheth et al 
20011 ) give a simple numerical fit for the impact 
of ellipticity on the linear collapse threshold, 6 ec : 



S ec Oi S sc [l + a(M) 



1.231 



(17) 



where S sc ~ 1.686 is the spherical collapse thresh- 
old. Since Sl oc t 2 / 3 in an Einstein-de Sitter uni- 
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verse, this implies that, for the same initial over- 
density, the collapse time of an elliptical region is 
longer than that of a spherical region by a factor 
of: 



t collapse, elliptical 



1.2313/2 



t 



collapse, spherical 



[1 + cr(M) 1 ' 23 ] 



(18) 



In other words, we need to extrapolate the linear 
theory predictions in Figs. farther beyond 

the point of intersection to actually hit gravita- 
tional collapse. Therefore, assuming that Equa- 
tion (|18l) is approximately independent of the dark 
energy model, it can be combined with the results 
of previous sections to show that the point of min- 
imum variance in 5l is shifted to smaller values of 

t/icollapse: i«e.. 



0.94 



collapse [1 + cr(M) 1 -' 



,3/2' 



(19) 



if we include the impact of ellipsoidal collapse. 
This result is shown in Fig. , and demonstrates 
how measuring the mass function of clusters at a 
given era may tell us about the entire history of 
linear growth, and not just a snapshot at the time 
of observation (as is in the traditional universal 
mass function hypothesis). 

Inspired by the fitting formula Eq. (|3|), Equa- 
tion (|19|) leads us to propose a new universal clus- 
ter mass function: 



f(a;Q) ~ g(a)e 



-h(a)/al 



, er* = oD 



0.94 x t a 



(1 



S)3/2 



(20) 

where the actual mass function is related to 
f(a; Q) through equation (TTJ, and t a is the time of 
observation at which D(t) is normalized to unity. 
In other words, Equation (|20]l suggests that the 
exponential cut-off in the cluster number counts 
at any time t is set by the linear density fluctu- 
ations cr» at an earlier time t*, which is set by 
Equation ([19)) . As suggested by several recent nu- 
merical studies (see Introduction) /(er; Q) depends 
on cosmology (denoted by Q), but we propose the 
functions g and h to be universal, while the de- 
pendence on cosmology (or dark energy models) 
only enters through D(t*). The explicit depen- 
dence of g and h on a at the time of observation is 
justified, as the value of a acts as a p roxy for the 



asphe ricity of the collapsing region (jSheth et al 
l200lh . 



Comparing to Eq. ([3]), we fix g and h as: 
i * . 1 . f 0.94 x t 

ffto = [UJ + !] > h (°) = CD [ (1 + (T 1.23)3/2 

n (21) 

where (A, a, 6, c) = (0.186, 1.47, 2.57, 1.19)0. 

Fig. ([5]) compares the evolution in the function 
f{o", Q) in N-body simulations, with our predic- 
tion from Eqs. (|20j|21l) . We see that while the 
result from spherical collapse, t* = 0.94 x t co u paae , 
underpredicts the evolution (top panel), Eq. (|20l) . 
which includes the ellipsoidal correction to the col- 
lapse time (Eq. [19]) , can successfully reproduce the 
evolution in /(cr; Q) (bottom panel). 

7. Conclusions and Future Prospects 

We have presented a study of non-linear grav- 
itational collapse in different models of dark en- 
ergy/cosmic acceleration. In particular, we crit- 
ically examined the correlation between the lin- 
ear growth of fluctuations and the emergence and 
statistics of collapsed objects (such as dark mat- 
ter haloes or galaxy clusters). First, we focused 
on the collapse of spherical overdensities, and dis- 
covered that they all have the same linear ovder- 
density (~ 1.50), at ~ 94% of the time of col- 
lapse/virialization, independent of the density or 
equation of state of dark energy. We then used a 
simple prescription to include the impact of el- 
lipsoidal collapse in our finding, and then used 
this result to propose a new universal mass func- 
tion for galaxy clusters/dark matter haloes. Our 
semi-analytic predictions are consistent with the 
observed evolution in mass function of haloes in 
N-body simulations. 

Future work will include the adaptation of this 
prescription to study models with dynamical equa- 
tions of state, such as quintessence or modified 
gravity models. A particularly challenging appli- 
cation woul d be to the gravitational aether /black 
hole model ( Prescod-Weinstein et al.l [20091) . Be- 
cause of the way dark energy is sourced in the 
gravitational aether model, there are unique tech- 
nical challenges associated with properly describ- 
ing structure formation in its presence. 



4 We recognize that these choices for g and h are not unique, 
but as Fig. (0 demonstrates, these are sufficient to fit 
current simulations 
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Fig. 6. — The relative change in /(<x) at high rcdshifts, 
compared to z = 0. The colo r ed re gions show the simu- 
lated results from lTinker et aL I d2008t> . Curves in top panel 
show our analytic prediction without ellipticity corrections, 
while curves in the bottom panel include the ellipticity cor- 
rections (Eq. I20|l . The solid, dotted, short-dashed, and 
long-dashed curves refer to z = 0, 0.5, 1.25, and 2.5 respec- 
tively. 



Finally, we recognize that a more systematic ap- 
proach to the question of universal mass functions 
should be possible, and given the level of theoret- 
ical and observational activities in this field, it is 
unlikely that the present work provides the last 
word on this subject. However, the novel (albeit 
trivial) observation of this work is that measur- 



ing the cluster mass function will teach us about 
the entire history of linear growth. This is in con- 
trast to many previous cosmological applications 
of cluster mass functions, which assume a one- 
to-one correspondence with the linear a(M) at 
the time of observation. Similar to a multi-level 
archeological excavation, dark matter haloes can 
be thought of as artifacts of linear growth. As Fig. 
(JSJ and Eq. ([20]) suggest, the number of smaller 
haloes (with larger er(M)) can teach us about the 
early evolution of linear growth, while the bigger 
haloes (with smaller a(M)) tell us about its more 
recent history. 

We thus speculate that this perspective can 
eventually lead to yet untapped information about 
the nature and history of cosmic acceleration, es- 
pecially as the releases of several large scale clus- 
ter surveys s uch as Atacama Cosm ology Tele- 
scope (ACT) dMenanteau et al.ll201Cl). So uth Pole 
Telescope (SPT) (lAndersson et alj |2010l). Planck 
(|Geisbusch k Hobsonl l2007t) . and Red Sequ ence 



Cluster Sequence 2 (RCS2) (|Yee et all 120071 ) are 
now looming on the horizon. 
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